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Abstract 

It is shown that a variational approach with fermionic squeezed states to many- 
fermion systems such as pairing model is one of useful methods beyond the usual 
Hartree-Fock-Bogoliubov approximation. A pairing-type quasi-spin squeezed state is 
constructed and adopted as a trial state in the variational method. By using this state, 
quantum fluctuations are taken into account. In pairing model, dynamical and static 
approaches are discussed in the context of comparing with the Hartree-Fock-Bogoliubov 
approximation about ground state energy. As a possible extension to the 0(4) model, 
prospective squeezed states are investigated and the usefulness is discussed. 
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1. Introduction 



As a possible extension of the time- dependent Hartree-Fock (TDHF) theory, a quasi- 
spin squeezed state was introduced instead of the Slater determinantal state as a possible 
extension of the trial state of variation. EM- 1 This state realizes the minimum uncertainty 
relation as well as the Slater determinantal state as a certain kind of coherent state. However, 
contrary with the coherent state approach, the quasi-spin squeezed state approach can take 
account of degrees of freedom for quantum fluctuations dynamically. Hereby, the squeezed 
state approach gives the better approximation than the coherent state approach in general. 

In the preceding work,B'II* applications of the quasi-spin squeezed state for many-fermion 
systems carried out in simple algebraic models. For example, the application to the Lip- 
kin model!-* where the particle-hole interaction is active shows clearly that the quasi-spin 
squeezed state gives the useful and powerful approximation, & because quantum fluctuations 
are contained properly. In this model, the quasi-spin squeezed state was constructed so as to 
take account of the particle-hole correlation. We call here it the particle-hole type squeezed 
state. This state reproduced the phase transition neatly, and in a certain limiting case 
for this state, the RPA equation was recaptured. 0*'i- ) From this viewpoint, the variational 
method with the quasi-spin squeezed state surpasses the RPA for many-fermion systems. 

This paper is devoted to the purpose of the application and promotion of the quasi-spin 
squeezed state to the models with pairing interaction. One of the present authors (Y.T.) 
with Yamamura and Kuriyama have given a way to construct the quasi-spin squeezed state 
for pairing model. B* Also, two of the present authors (H.A. and Y.T.) have given a possible 
treatment of the degree of freedom for the quantum fluctuations in the Lipkin model. 
In this paper, we extend the treatment previously given & to the dynamical treatment in 
the pairing model. In the pairing model, we investigate the ground state energy with two 
different manners using the quasi-spin squeezed state. One is that the quasi-spin squeezed 
state is used in the framework of the time-dependent variational principle, and the system is 
considered dynamically. The other is that the energy minimum is sought directly introducing 
the chemical potential for the particle number conservation. 

Further, we extend the quasi-spin squeezed state for the pairing model to the 0(4) model. 
Recently, to investigate the shape coexistence phenomena of nucleus, the simplest 0(4) 
algebraic model with the pairing and the quadrapole interactions has been reanalyzed beyond 
the Hartree-Fock-Bogoliubov and the random phase approximations. As an extension 
of the paring type quasi-spin squeezed state, we attempt to construct possible quasi-spin 
squeezed states for the 0(4) model. 

This paper is organized as follows. In the next section, the exact treatment and the 
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Hartree-Fock approximation of pairing model are recapitulated containing notation. 

In §3, the time-dependent variational approach with the quasi-spin squeezed state for the 
pairing model is formulated. Also, it is shown that the ground state energy is well reproduced 
analytically in this dynamical approach. Against the former section, in §4, a static treatment, 
that is, the variation for the Hamiltonian with the chemical potential for the particle number 
conservation is carried out, and the ground state energy is estimated numerically in the 
pairing model. In §5, the extensions of the pairing-type squeezed state to the one for the 0(4) 
model are given. The energy expectation value for the ground state is calculated numerically 
by using two candidates of the squeezed states for the 0(4) model. The expectation values 
for the various operators with respect to two squeezed states are summarized in Appendix 
A and B, respectively. The last section is devoted to a summary. 

§2. Recapitulation of exact solution and coherent state approximation for 

pairing model with single energy level 

In this section, the exact treatment of an exactly solvable quantum many-fermion model, 
which is called the pairing model, is reviewed for later convenience. Also, the coherent state 
approximation, which corresponds to the BCS approximation, is recapitulated. 

2.1. Exact solution for the pairing model 

We investigate a simple many-fermion system in which there exists N identical fermions 
in a single spherical orbit with pairing interaction. The single particle state is specified by a 
set of quantum number (J, m), where j and m represent the magnitude of angular momentum 
of the single particle state and its projection to the z-axis, respectively. Thus, let us start 
with the following Hamiltonian : 

H = eJ2clc m -^J2(-) i ~ m ^J-mj:(-) j - m C- m &m , (2-1) 

m m m 

where e and G represent the single particle energy and the force strength, respectively. The 
operators c m and & m are the fermion annihilation and creation operators with the quantum 
number m, which obey the anti-commutation relations : 

{ Cm > ^m' } = 3mm' ) { C m , C m i } = { , C m i } = . (2'2) 

We introduce the following operators : 

s + = \Y.(-r m t ] J~m, 5_ = V(-r^ m , s = Wcic m -n) , (2-3) 

" m m m 
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where ft represents the half of the degeneracy :ft = j + 1/2. These operators compose the 
sw(2)-algebra : 

[S + ,S_] = 2S , [S Q ,S ± ] = ±S ± . (2-4) 



Thus, these operators are called the quasi-spin operators. eD'ESP Then, the Hamiltonian (|2-1 ) 
can be rewritten in terms of the quasi-spin operators as 

H = 2e(S + Sj) - GS+S- 

= eN- GS+S- , (2-5) 
Si = ft/2 (= (j + l/2)/2) , 

where N represents the number operator : 

# = E4*»- (2-6) 



ni 



As is well known, the eigenstates and eigenvalues for this Hamiltonian are easily obtained. 
From [So, H] = 0, the eigenstates are given in terms of the eigenstates of So : 

S \S, Sq) = S(S + i)\S, Sq) , Sq\S, So) = Sq\S, So) , 



S±\S, So) = y/(S =F S )(S ±S + 1)\S, So ± 1) , (2-7) 

where S = Sq + (S + S^ + S-.S+)/2. Thus, we obtain the eigenvalue equation 

H\S, So) = [2e(S + S,) - G(S + S )(S -S + 1)]\S, S ) = E V \S, S ) . (2-8) 

Here, from the relation So = N/2 — Sj, the eigenvalue So can be written in terms of particle 
number iV as So = N/2 — Sj. Further, if we introduce the variable u as 

S = Sj - u/2 , v = 0, 1, ■ • • , 2Sj (= ft) , (2-9) 

the energy eigenvalue E u is given by 

E u = eN - ^GN(1 - u/N)(2ft - N - v + 2) . (2-10) 

Thus, the ground state energy can be obtained by setting v = as 
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2.2. Coherent state approximation 

Next, we review the coherent state approach to this pairing model, which is identical 
with the BCS approximation to the pairing model consisting of the single energy level. 
The sw(2)-coherent state is given as 

10(a)) = U(f)\0) , 

U(f) = ex V (fS+ - /*&-) , 

|0> = \S = Sj, S = -Sj) , (5_|0) = 0) . (2-12) 
We impose the canonicity condition : 

W«)I^W«)> = \C , = ~i • (2-13) 

A possible solution of the above canonicity condition is presented as 

f • I £ I I c\ ~C1 f 



i = \l*SiTZ sin \f\ . C = sin |/| . (2-14) 



Then, the expectation values are obtained as 



(<f ) (a)\S + \<J ) (a))=epS 3 -\ZW 
(<f>(a)\S-\<f>(a)) = y/2S j -\£\*£, 

(<j ) (a)\S \<j ) (a)) = \^-S j , (2-15) 
Ma)|£ + S-W«)> = 2^-lCl 2 " ( X " ^r) 1*1* > (2' 16 ) 

Thus, the expectation value of Hamiltonian H and number operator N are given as 

(<P{a)\H\<P{a)) = e ■ 2|£| 2 - G faltf - (l - |£| 4 ) = E , 

(0(«)|iV|0(a)) = 2|e| 2 = iV. (2-17) 

If total particle number conserves, that is, N =constant, then, the energy expectation value 
E is obtained function of N : 

1 / N N \ 
E = eN-- 4 GNn 2-- + s5 . (2-18) 



It is interesting to compare the exact ground state energy ( 2T1|) and the BCS approx- 



imated energy ( 2-18| ). The last term of ( j2-18 ) is different from that of the exact result in 



( |2 • 1 1| ) . Let us assume that the particle number N and the half of the degeneracy Q are the 
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same order of magnitude. If Q (or N) is large, both the last terms of the exact eigenvalue 
( 2- 11|) and the approximate ground state energy ( 2-18|) can be neglected. Thus, the coher- 
ent state approximation presents good result for the ground state energy. This situation is 
similar to large N limit which occurs in the several field of physics. 

In the next section, we try to reproduce the last term in a certain approximate approach, 
namely, the time-dependent variational approach with quasi-spin squeezed state. 



§3. Dynamical approach of quasi-spin squeezed state for pairing model 

In this section it is shown that the ground state energy of a many-fermion system with 
pairing interaction can be well approximated analytically by using the time-dependent vari- 
ational approach with a quasi-spin squeezed state. 

3.1. Quasi-spin squeezed state 

In this subsection, the quasi-spin squeezed state is introduced following to Ref.p. First, 
we introduce the following operators : 

i f = ^=, ^=4=r> N = 2(S + Sj), (3-1) 
2Sj 



where N is identical with the number operator ( |2-6|) . Then, the commutation relations can 
be expressed as 

[A, it] = i _ A , [n, A] = -2A , [N, it] = 2 it . (3-2) 

Using the boson-like operator A\ the sw(2)-coherent state in (2-12) can be recast into 

10(a)) = I = exp(ait)| ) ) 

<P(a*a) = 1 1 + |-J , (3-3) 

where a is related to / in (2-12) as a = j2Sj ■ (f /\f\) ■ tan |/|. The state 10(a)) in (3-3) is 
a vacuum state for the Bogoliubov transformed operator a m : 

a m = Uc m - V(-y- m cl m , a m \<f)(a)) = . (34) 

The coefficients U and V are given as 

U= , 1 V= " . 1 U 2 + \V\ 2 = 1. (3-5) 
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Of course, a m and are fermion annihilation and creation operators and the anti-commutation 
relations are satisfied. By using the above Bogoliubov transformed operators, we introduce 
the following operators : 

^ t = ^E(-r m ^ f - m , 5 = -i=5:(-r m a_ m a m , M = ^ala m . (3-6) 

\/OOj m \/(iOj m m 

Then, the state \<f>((x)) satisfies 

B\(j){a)) = M\(j)(a)) = . (3-7) 
Further, the commutation relations are as follows : 

[B, &} = 1 - ^- , [M, B] = -IB , [M, fit] = 2 fit . (3 . 8 ) 

The quasi-spin squeezed state can be constructed on the sw(2)-coherent state |0(a)) by using 
the above boson-like operator B^ as is similar to the ordinary boson squeezed state : 

^(a, (3)) = 1 exp (\PBV) |0(a)) , 

We call the state \i/j(a, (3)) the quasi-spin squeezed state. 

3.2. Canonicity conditions and expectation values 

In this subsection, we calculate the expectation values for various operators. The ex- 
pectation values can be expressed in terms of the canonical variables which are introduced 
through the canonicity conditions. The same results derived in this subsection are originally 
given in the Lipkin model in Ref.0) at the first time and these results were used in Ref.|6]) 



to analyze the effects of quantum fluctuations in the Lipkin model. The Lipkin model is 
a simple many-fermion model with two energy levels in which the particle-hole interaction 
is active.!^ This model has the same algebraic structure as the pairing model, that is, the 
Hamiltonian can be expressed by the quasi-spin sw(2)-generators. Thus, in our case, the 



same calculation as that given in Ref.|T2|) are carried out. 



For the quasi-spin squeezed state \ip(a,f3)) in (3-9), the following expression is useful : 

+ i 1 ~ —¥tm ) ' 1 + <**«/(2Sj) 2 {a d ' a ~ ^ } ' 
nu) S *f£> , (3-10) 



7 



where d z = d/dz. The canonicity conditions are imposed in order to introduce the sets of 
canonical variables (X, X*) and (Y,Y*) as follows : 



(tlj(a,(3)\d x \^(a,(3)) = -X* 



(i;(a,P)\dy^(a,p)) 



y* 



(1>(a,P)\d x -\iJ>(a,P)) = -~X 
(ij>(a,l3)\dy.\iP(a,P)) = ~Y . 



(3-11) 



Possible solutions for X and Y are obtained as 

X 



X* 



a 



P 



1 



X*X 4Y*Y 



a 



2Si 



25, 



x *x 4Y*Y 



2Sj 



25, 



Y 



P* 



Y* 



/ K(Y*Y) ^K(Y*Y) 
where K(Y*Y) is introduced and satisfies the relation 

K(Y*Y)W(Y*Y/K) = V'(Y*Y/K) . 



(3-12) 



(3-13) 



The expectation values for B, B\ M and the products of these operators are easily obtained 
and are expressed in terms of the canonical variables as follows : 



(^{a,p) 
(ip(a,f3) 

(ip(a,/3) 

(ip(a,P) 
(ip(a,p) 



B\1>(a,0)) = (,f>(a,P)\&Ma,P))=O 
M\^{a,P)) = 4Y*Y , 



B 2 ty(a,P)) = 2YJK(y*Y) 



B^(a,P)) = 2Y*^/K(Y*Y) , 

&B\i>(a,P)) = 2(1 - 1/{2Sj))Y*Y - ^(Y*Y) 2 ■ L(Y*Y) 

BB^(a,P)) = (il>(a,P)\&B\il>(a,P)) + (1 - 2Y*Y/S j ) , 
M 2 \ijj(a, P)) = 16Y*Y(1 + Y*Y ■ L(Y*Y)) , 



(3-14a) 
(3-14b) 

(3- 14c) 
(3-14d) 



where L(Y*Y) is defined and satisfies 



K(Y*Yy ■ L(Y*Y) 



V"(P*P) 

*Wp) 



(3-14e) 



By using the relations between the original variables a and P and the canonical variables 
X and Y, the coefficients of the Bogoliubov transformation (|3-5|), U and V, are expressed as 



U 



1 - 



X*X _ 4Y*Y 

2 jS" ri 2 S A 



1 - 



4Y*Y 
25,- 



V 



X 



V 1 4 25, y 



(3-15) 
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Then, the operators A, A> and N, which are related to the quasi-spin operators S-, S + and 
So, respectively, in ( |3-1| ), can be expressed as 

M 



A = J2SiUV 1 - 



2S. 



- V 2 B ] + U 2 B , 



A ] = JlSjUV* ( 1 - ^-J + U 2 & - V* 2 B , 

N = ASjV*V U - + J2S~U{V& + V*B) + M . (3-16) 

Thus, the expectation values for A, A\ N and the products of these operators are easily 
obtained and are expressed in terms of the canonical variables as follows : 



y>{<x,p)\A\4(a,P)) = X*J\ 
(^(a,(3)\A^ma,(3)) = 



X*X AY*Y 



2 Sj 2 S j 



s X*X 4Y*Y 
2 Sj 



(ip(a, 0)\N\ip(a, (3)} = 2X*X + 4Y*Y , 
(iP(a,(3)\A 2 \iP(a,P)) = ^{a,(3)\A\^{a,(3)) 2 

U 2 V 2 



(3-17a) 



2Sj - AY*Y 



2S< 



^{a,(3)\A^{a,(3)) = {^{a, (3)\A^{a, (3)) 2 

U 2 V* 2 



(^(a,f3)\M 2 \^(a,(3)) - (^(a, (3)\M\^(a, f3)) 2 

+V 4 (^(a, (3)\B^ 2 ma, (3)} + U 4 (^(a, (3)\B 2 \^(a, /?)) 

-2U 2 V 2 ^(a,{3)\&B\^a,{3)) , 

1 



2Si - AY*Y 



2Si 



(^(a,(3)\M 2 \ij(a,/3)) - (*P(a, /3)\M\if>(a, f3)) 2 _ 

+U\i,{a, {3)\B^(a, (3)) + V*\^{a, (3)\B 2 \^(a, (3)) 
-2U 2 V* 2 (tfj{a,(3)\B^B\tfj{a,p)) , (3- 17b) 



(iP(a,P)\AiA\i/>(a,P)) = P)\A^(a, P))(iP(a, f3)\A\iP(a, (3)) 

{X*Xf 



2Sj(2Sj-AY*Y) 
U 2 V*V r 



2S 



(tfj(a,(3)\M 2 ma,(3)) - (</>(«, /3)\M\4>(a, (3)) 2 

U 2 V 2 (i;(a,P)\B^(a,P)) - U*V* 2 {^{a, (3)\P\^{a, (3)) 
{l-2U 2 V*V)^{a,(3)^B\^{a,(3)) , 



(^{a, (3)\AA^(a, (3)) = (^{a, (3)\A^A\^(a, (3))+il 



2X*X 4Y*Y 



2Sj 2S 3 



(3-17c) 
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From (|3-1|), the expectation values of quasi-spin operators are derived from (3-17a) : 

(ip(a,/3)\S+\ip(a,/3)) = X*^2Sj - X*X - AY*Y , 
(ip(a,P)\S-\ip(a,P)) = ^25,- - X*X - 4Y*Y X , 
(tj)(ai, P)\So\tJ;(a, 0)) = X*X + 2Y*Y - Sj . (3-18) 

By comparing with (2-15), the variable 2\Y\ 2 represents the quantum fluctuations. 

3.3. Time- dependent variational approach with quasi-spin squeezed state for pairing model 

The model Hamiltonian (2-5) can be expressed in terms of the fermion number operator 
N and the boson-like operators A and A) as 

H = eN-2S i GA^A. (3-19) 



Thus, the expectation value of this Hamiltonian is easily obtained by using Q3-17 ). We denote 



it as H sq 



H sq =(iP(a,f3)\Hma,f3)) . (3-20) 



The dynamics of this system can be investigated approximately by deriving the time- 
dependence of the canonical variables (X, X*) and (Y, Y*). The expectation value for the 
time-derivative is calculated as 

<V(a, P)\idMa, (3)) = % -{X*X - X*X + Y*Y - Y*Y) , (3-21) 

where X represents the time-derivative of X and so on. The time-dependence of X, X*, Y 
and Y* is derived from the time-dependent variational principle : 



5 J (ip(a,P)\id t - H\ip(a,p))dt = . (3-22) 

Hereafter, we assume that \Y | 2 <C 1 because \Y | 2 means the quantum fluctuations. Then, 
K(Y*Y) and L{Y*Y) defined in ( |3-13|) and ( |3-14e|) , respectively, can be evaluated by the 



expansion with respect to \Y\ . As a result, we obtain 

^=K 1 -i) + ( 1 -i + 4F) 5 ^ + 

i 



9<?- / 
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Further, we introduce the action-angle variables instead of (X,X*) and (Y,Y*) as 



'n x e 



Y = ^e~ idY , Y* = ^e l0Y . 



(3-24) 



Then, the expectation values for the Hamiltonian, the time-derivative and the number op- 
erator can be expressed as 



2S jnx - n x + ^ 



H sq = e(2n x + An Y ) - G 

-2^1 - ^- (l - n x ^cos{29 x - 6 Y ) 



10 2 2 8 



+2 ( 2S j - 1 - An x + — n x + —n x - 



2S, 



(i/j(a,/3)\id t \iP(a,(3)) = (n x 9 x + n Y 9 Y ) , 
iV = (i/}(a, (3)\N\if)(a, (3)) = 2n x + An Y . 



(3-25a) 
(3-25b) 
(3-25c) 



From the time-dependent variational principle ( |3-22| ), the following equations of motion 
are derived under the above-mentioned approximation : 



e 



dH. 



x 



sq 



dn x 



2e-G-2 



Sj-n x - 



n x 



2S, 



V2^1 - X^i _ M^cos(2e x - Or) 

(3-26a) 



+4 T 1 + IS" + 2S™ X ~ Sj nx ) ny 



dK 



n x 



Sq 



-G • 4^1-7-^(1 - ^^n xv ^7sin(2^ - 6 Y ) , (3-26b) 



sq 



dm 



Ae-G 



■V2jl-J^( 1 - ^ N )-^cos(2^ - 



% \_nx_ 
2Sj\~ 2Sj)^w 

( 5 12 

+2 -1 + 2S j - An x + — + — n\ - —rn 2 



Q ' Q -X Q2'°X 



OK 



n Y 



sq 



86y 



G • 2^1-^1 - ^n x ^sm{29 x - 6 Y ) . 



(3-26c) 
(3-26d) 



It is found from (3-26b) and (3-26d) that the total fermion number N in (3-25c) is conserved, 
that is, 

N = 2n x + Ah Y = . (3-27) 
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3.4. Dynamical approach to the ground state energy 

It should be noted here that nx < N/2 from (3-25c) and N < Q/2 = Sj. Thus, the 
inequality 1 — n x /2Sj > is obtained. From the approximated energy expectation value 
(3-25a) for G > 0, the energy minimum is then obtained in the case cos(2#x — 6y) = — 1, 
namely, 

Q Y = 26 x + vr . (3-28) 

Since the energy is minimal in the ground state, the relation ( |3-28 ) should be satisfied at any 
time. In order to assure the above-mentioned situation, the following consistency condition 
should be obeyed : 

9 Y = 29 x . (3-29) 

Thus, from the equations of motion (3-26a) and (3-26c), under the approximation of small 
the consistency condition ( |3-29| ) gives the following expression of ny in the lowest order 
approximation of ny '■ 



1 - -i- 

25, 

ny ~ — T=7 V 

2^(1 - f 



i -i 

1 1 

1 + 



2r* (!-*£)(! 



(3-30) 



Thus, by substituting (|3-28|) and y/ny in (3-30) under the lowest order approximation of ny 
into the energy expectation value (3-25a), and by performing the approximation of large N 
or large Q{— 2Sj) approximation, we obtain the ground state energy as 

H sq = eN- hjNQ (2 - ^ + ^ + 0(l/JV/2, l/N 2 )^ . (3-31) 

This result reproduces the exact energy eigenvalue ( [2-18 ) by neglecting the higher order 
term of l/NQ, l/i7 2 and 1/N 2 for large iV and Q limit. Thus, the quasi-spin squeezed state 
presents a good approximation in the time-dependent variational approach to the pairing 
model. In this approach, the existence of the rotational motion in the phase space consist- 
ing of (nx, Ox'-, n Y) 9y) plays the important role. The angle variables for rotational motion, 



6x and 6y, are consistently changed in ( 3-29 ). This consistency condition is essential to 
reproduce the exact energy for the ground state under the large iV and i? limit. The approx- 
imation corresponds to so-called large iV approximation. In general, it is known that the 
large N expansion at zero temperature corresponds to h expansion. In this sense, the time- 
dependent variational approach with the quasi-spin squeezed state gives the approximation 
including the higher order quantum fluctuations than h if any expansion is not applied. 
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§4. Numerical estimation of ground state energy using quasi-spin 

squeezed state 

In the previous section, it has been shown that the exact ground state energy for the 
pairing model can be well reproduced in the large N and Q approximation in the time- 
dependent variational approach with the quasi-spin squeezed state. In that treatment, the 
rotational motion, which is originated from the use of the number violating state, plays the 
essential role to reproduce the ground state energy. In this section, we evaluate the expec- 
tation value for the ground state energy numerically by using the quasi-spin squeezed state 
without the expansion of N and Q. Instead the time-dependent variational approach with 
the consistency condition developed in the previous section, the fermion umber conservation 
is guaranteed by introducing the chemical potential. 

We impose the minimization condition for the following quantities : 

6(H') = 8{H - fiN) 

= ((e - fi)N - 2S 3 GA^A) = , (4-1) 

where /1 represents the chemical potential and (• • •) denotes the expectation value with 
respect to the quasi-spin squeezed state. The expectation values for N and A^A have been 
already given in (|3-17|) . The variation can be carried out the variational parameters (a, a*) 
and If we put (3 = (3* = 0, the state is reduced to the sw(2)-coherent state. 

In Fig.l, the ground state energy with the unit e is depicted in the case N = f2 = 8. The 
horizontal axis represents the force strength G of the pairing interaction. The dotted curve, 
dot-dashed and solid curves represent the exact energy eigenvalue, the expectation value of 
the Hamiltonian with respect to the coherent state and the quasi-spin squeezed, respectively. 
The result of the quasi-spin squeezed state well reproduces the exact eigenvalue for the wide 
range of G, comparing with the result by the coherent state. 

In Fig. 2, the energy is depicted in the case e = 1.0 and G — 1. The horizontal axis 
represents the particle number iV with Q = N. The result obtained by using the quasi-spin 
squeezed state is almost same as the exact eigenvalue. These figures show that the squeezed 
state approach presents a good approximation. 

§5. Extension to 0(4) model with pairing plus quadrapole interactions 

In the previous sections, §§3 and 4, the sw(2)-algebraic model with the pairing interaction 
has been investigated by using the quasi-spin squeezed state. It has been shown that the 
ground state energy has been well reproduced compared with the usual sw(2)-coherent state. 



13 




0.5 1 1.5 2 



G/e 

Fig. 1. Energy expectation values with respect to the coherent state (dot-dashed curve) and the 
squeezed state (solid curve) are depicted together with the exact eigenvalues (dotted curve) in 
the case N = ft = 8. The horizontal axis represents G with the unit e. 
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Coherent 
Squeezed 
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Fig. 2. Energy expectation values with respect to the coherent state (dots) and the squeezed state 
(diamonds) are depicted together with the exact eigenvalues (crosses) in the case e = 1.0 and 
G = 1.0. The horizontal axis represents N. 



In this section, we try to extend our squeezed state approach to the 0(4)-algebraic model 
with both the pairing and the quadrapole interactions in the many-fermion system such as 
nucleus. 



5.1. 0(4) model with pairing and quadrapole interactions 

Let us start with the single-j shell model, where j represents the angular momentum 
quantum number. Thus, the degeneracy 2i? is 2i? = 2j + 1. The pairing and the quadrapole 
interactions are active in this model. The Hamiltonian can be expressed as© 



(5-i; 
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Here, we define the following operators in terms of the fermion annihilation and creation 
operators {c m , c^} as 



c m c m , 



^=E44, P=± 

m>0 m>0 
3 ^3 

Q = (7 rnC rn C m , N = C m C m , 

m=—j m=—j 



3 3 

it ^ 

m>0 m>0 



P f = £ ffm4c~ , P = E ^mC~C m , (5-2) 



where 



( 1 , for m < 
m y ' m m \ _i for m > /2/2 



and J? represents the half of the degeneracy. In (5-2), we also define and P for the later 
convenience. 

The Hamiltonian (|5-1|) has the 0(4)-algebraic structure. We can construct two su(2)- 
generators from the operators in (5-2) : 

$5. = (pt + j*) = e 4^ = 0^, 

0<m<Q/2 

Sj = (iV + Q-^)/4 = (2 E 4^-^1/4, (54a) 
= = E 44 = (5?)t, 

i?/2<m<r2 

5 ° = (iV-g-/2)/4= [2 E 4^-^]/4, (54b) 

\ i?/2<|m|<« / 

where these operators satisfy the following commutation relations : 

[5+7 ] = 2S"q , [iSq) = i>S± j 

[5°, 5°] = 24 , [s°,s2] = ±s2, 

[Si,o, S2, ] = • (5-5) 

Thus, the sets {£±0} an d {S±o} gi ye two sets of independent sw(2)-generators. By express- 
ing the operators P\ P and N inversely in terms of the above two sw(2)-generators, the 
Hamiltonian in ( |5-1| ) can be expressed in terms of two independent sw(2)-generators. Thus, 
this model given by the Hamiltonian in ( |5T| ) is called 0(4) (~ su(2) x su(2)) model. This 
model can be solved exactly because the model space is spanned by two quasi-spin su{2) 
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states and the diagonalization is easily performed. Thus, the validity of an approximation 
can be checked. 

An approach by the 0(4)-coherent state corresponds to the Hartree-Fock approximation. 
The 0(4)-coherent state can be constructed by the direct product of the two su(2)-coherent 
state as 

|c) (4) = \oci) <g> Wo) = \a I} a ) , 

1 / §i \ 1 f S° \ 

Wl) = ~7^= eX P a I /cTo~ 1°)/ > Wo) = ^= ex P a O 7=j= |0)o , 



Sj = Q/A , S = n/A . (5-6) 

Thus, the variational parameters are aj, a}, ao and a* in 0(4)-coherent state. 

For the number conservation, we introduce the chemical potential. Then, the variation 
with respect to a/, a}, ao and a* Q is carried out as 

5(ai,ao\Ho(4) — fiN\ai,ao) = ■ (5'7) 

Figure 3 shows the exact ground state energy eigenvalues (dotted curve) and the energy 
expectation value calculated by the 0(4)-coherent state (dashed curve). The model parame- 
ters G and x are parameterized as G = cos 9 and 2\ = sin 9, respectively,^ and iV = Q = 8 
and e = 0. It is found that, if the quadrapole interaction is dominant, that is, 9 is large, the 
coherent state approximation presents good results for the ground state energy. However, 
the pairing interaction is dominant, the coherent state approximation is rather bad. There is 
a room to devise the approximate state in the 0(4) model with the pairing and quadrapole 
interactions such as nucleus. 

5.2. Direct product of two quasi-spin squeezed states 

As is similar to the su(2)-squeezed state given in §3, let us construct the squeezed state 
for the 0(4) model. First, let us introduce the Bogoliubov transformed fermion annihilation 
and creation operators (a m , b m ; a^, fetj : 

a m = Uic m - V/cL , for \m\ < Q/2 , 

b m = U c m - V ct , for Q/2 < \m\ < fi , (5-8) 
where the coefficients of the Bogoliubov transformation are given as 



1 OtT 1 1 : 

Ui= Vj= _L_ uf+ v 

'I + a* I a I /(2S I ) VZSj i + a * I a I /(2S I ) 
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Uo= , V = 1° . Ul + \V\l = l 

'l + a a /{2S ) V2So Jl + a a /(2S ) 



(5-9) 

Then, the boson-like operators are introduced like ( |3-6|) : 
1 1 

B ^ = ~7W ' B I = -/f&= E U~a m , Mj= E h \rflm , 

0<m<i?/2 V^^I <m <n/2 \m\<Q/2 

b ] = -i E tiJl . = E hbm , M = e • 

(5-10) 

Then, the following commutation relations are satisfied : 

[B T , B\] = l-^r, [Mi, Bi] = -2B T , [Mi, B\] = 2B\ , 

[B , B Q ] = 1 - ^ , [Mo, S ] = -25 , [Mo, B ] = 2B . (5-11) 

The quasi-spin squeezed state \ip) for 0(4) model may be constructed by the direct 
product of the two quasi-spin squeezed states, \ip{ai, Pi)) and \i[)(ao, Po)), which are defined 
as is similar to (3-9) : 

\i/>) = \i){ai,(3i))® \i){a ,Po)) , 

1 /I - A t 2 



\ip(ai,Pi)) = 7= exp -^/-Bj |a/) , 



1 /I - A t 2 



|^(ao,/3o)> = ^====exp ( -/3 5 ) |a G ) , (5-12) 



*Wo) = i + e i ^L# : (i^i 2 ) fc n f 1 ; 



- (2*)" ' ^ V 25 G 

Since the two sw(2)-generators are expressed in terms of the above introduced boson-like 
operators in (5-10) as 

Si = 2Si (l - EW + ^i{U 2 jB\ - v; 2 B!) = (sly , 
S 1 , = -Si(U] - \Vi\ 2 ) (l - ^) + \fiSiUi(ViBi + V?B 
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Fig. 3. Energy expectation values with respect to the coherent state (dashed curve), the direct 
product of two quasi-spin squeezed state (dot-dashed curve) and the fermionic squeezed state 
(solid curve) are depicted together with the exact eigenvalues (dotted curve) in the case N = 
Q = 8. The horizontal axis represents 9 which parameterize the model parameters G = cos 9 
and 2x = sin 9. 



S° = 2S [l 



^ U VS + j2S~oiU 2 Bl - VS 2 B Q ) = (S°y 



2S, 



o , 



4° = -S {U 2 o - \V C 



o\ 



1 - 



]f ° 1 +J2S~ Uo{V Bo + V*Bl) 



2S. 



o , 



(5-13) 



Also, the pairing operator P, the quadrapole operator Q and the number operator N can 
be expressed in terms of the above two sets of the sw(2)-generators as 



p = Si + S° , Q = 2(Sl- 4°) , A> = 2{Sl + S °) + n . 



(5-14) 



Thus, the expectation value of the Hamiltonian ( |5T| ) can be expressed by the sw(2)-generators. 

The variation are carried out with respect to the eight variational parameters aj t o, &* I O , 
Pi t o and Pj as 

6{il>\H O (4)-iiN\il>) = 0. (5-15) 

The expectation values are summarized in Appendix A. The energy expectation value for 
this model Hamiltonian with respect to the direct product of the two quasi-spin squeezed 
state are depicted in Fig. 3 (dot-dashed curve) compared with the exact eigenvalues (dotted 
curve) and the expectation values derived by using the 0(4)-coherent state (dashed curve) 
in the case N — f2 — 8. The horizontal axis denotes 9 where we parameterize the force 
strength G and x as G = cos# and 2\ = sin#. It is found that the direct product of the 
two sw(2)-spin squeezed states does not present good results for the ground state energy, 
especially in the region where the pairing interaction is dominant. The reason is as follows : 
If x — 0, the 0(4) model is reduced to the pairing model investigated in §§2~ 4. However, 
the direct product of the two quasi-spin squeezed state is not reduced to the state in (3-9) 
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because of the absence of the cross term B}B , even if the parameters are set up as aj = ao 
and Pi = Po- Thus, in this case, the direct product does not include appropriate pairing 
correlations fully, and thus does not give a suitable squeezed state for many-fermion systems. 

5.3. Fermionic squeezed state 

In the previous subsection, the squeezed state for the 0(4) model has been constructed 
by the direct product of the two quasi-spin squeezed state. The extension of the sii(2)-quasi- 
spin squeezed state to that for the 0(4) model may be natural from the viewpoint of the 
algebraic structure. However, the ground state energy is not reproduced so well. In this 
subsection, the trial state can be devised, which we call a fermionic squeezed state. 

We define another squeezed state \s) as 

c) (4) , (5-16) 

= [(2n)!] 2 * OToTT! tt U _ M 2n rf 1 (\ - -±- 
h W h (2n-r)!r! \k\ 2Si) £ I 2S Q 

The state \s) can be reduced to the quasi-spin squeezed state if the conditions ai = ao 
and Pi = Po are imposed. Thus, the pairing correlations are fully taken into account. The 
necessary expectation values for calculating the expectation value of the Hamiltonian are 
summarized in Appendix B. Thus, we can derive the variational equations in Eq.( |5T5| ). By 
solving the variational equations, we can obtain eight parameters (aj,o, oiq o-> Pi,o, Pi o)i an d 
then the energy expectation value is evaluated. In Fig. 3, the energy expectation values 
calculated by using the state \s) are depicted (solid curve) compared with the energies of 
exact approach, coherent state approximation and the approximation by using the direct 
product of two quasi-spin squeezed state in (5-12). It is shown that our squeezed state 
approach with the state \s) presents good results. In Fig. 4, variables \Pi\ and \Po\, which 
represent the quantum fluctuations and by which the particle-particle correlations are taken 
into account, are depicted in the case of the fermionic squeezed state (a) and the direct 
product of two quasi-spin squeezed state (b), respectively. In both cases, \Pi\ (solid curves) 
and \Po\ (dashed curves) are almost the same. In all regions, the values \Pi\ and \Po\ in 
the fermionic squeezed state approach are larger than those in the direct product of two 
quasi-spin squeezed state approach. Especially, in the fermionic squeezed state approach, 
it is found that the values are not negligible in the region where the pairing correlation 
is dominant. Thus, the pairing correlation is taken into account in this state. However, 
the values of \Pi\ and \Po\ are rather small in the region where the quadrapole-quadrapole 
correlation is dominant. It may be concluded that the 0(4)-coherent state is rather good 
state for describing the system in which the quadrapole correlation is rather strong. 



\s) 



cxp 



PiB\ + P B ] 




19 




7T/4 7T/2 

13 = 1.0 a -o.o 

2 X = 0.0 2X = 1.0 



Fig. 4. The variables (solid curves) and \(3o\ (dashed curves) are depicted in the fermionic 
squeezed state approach (a) and the approach with the direct product of two quasi-spin squeezed 
state (b), respectively, in the case N = Q = 8. The horizontal axis represents 9 which param- 
eterize the model parameters G = cos 9 and 2\ = sin#. 
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Fig. 5. Energy expectation values with respect to the coherent state (dashed curve), the direct 
product of two quasi-spin squeezed state (dot-dashed curve) and the fermionic squeezed state 
(solid curve) are depicted together with the exact eigenvalues (dotted curve) in the case N = 
fi = 8 and x = 0-2- The horizontal axis represents G. 

In Figs. 5 and 6, the energy expectation values with respect to the coherent state 
(dashed curve), the direct product of two quasi-spin squeezed state (dot-dashed curve) and 
the fermionic squeezed state (solid curve) are also depicted together with the exact eigen- 
values (dotted curve) in the case N = fl = 8 with \ — 0-2 (Fig.5) and G = 1.0 (Fig. 6), 
respectively, in which the horizontal axes represent G (Fig.5) and x (Fig-6). It is found 
that the fermionic squeezed state approach presents rather good results in all the parameter 
regions with G and %. 
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Fig. 6. Energy expectation values with respect to the coherent state (dashed curve), the direct 
product of two quasi-spin squeezed state (dot-dashed curve) and the fermionic squeezed state 
(solid curve) are depicted together with the exact eigenvalues (dotted curve) in the case TV" = 
Q = 8 and G = 1. The horizontal axis represents \- 



5.4. k-squeezed state 

In this subsection, we give one comment. Let us consider the su(2) pairing model dis- 
^4. We can define the fc-squeezed state as 
1 / w 



cussed in 52 



\s)k 



Tfe exp (^-1^) |0(«)> 



(k > 3, (integer)) 



[ |? (kn) 



fen— 1 



4: n (n\] 



\w\ 



1 - 



V 



2Si 



2 n 

n=0 V-J p=l \ "^3 . 

In this state, the expectation values are calculated easily as 

k (s\B\s) k = k (s\B^\s) k = k (s\B 2 \s) k = fc (5|-B +2 |s)fc = 

P 



k (s\M\s) k = — J2 ( 2k 
1 1 



(n\y 4" U { 2S jy 



n 



k n =o 
[n/k] 



k (s\M 2 \s) k = — J2(2kn) 
1 1 



2 {kn)\ \w 



2n kn—1 



k n=0 



n! 



4" 



P =i 



1_ [ ^ ] {kn + l)\\w\ 2n ^ ( 



k ( s \BB^\s) k =-J2 

1 k n =0 

k (s\&B\s) k = k (s\BBi\s)k- ^1 



:(s\M\ S )k 

2Si 



(5-17) 



(5-18) 



and the other combinations of the products of two operators give no finite values. In order 
to calculate the ground state energy, we impose the condition of minimization of the energy 
expectation value. As a result, it is concluded that \w\ = is a solution for minimization. 
Thus, the fc-squeezed state is reduced to the su(2)-coherent state in this pairing model. 
Therefore, the quasi-spin squeezed state gives a good state in order to take into account the 
pairing correlation. 
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6. summary 



We have investigated a validity of the quasi-spin squeezed states as trial states in the 
variational method in the many-fermion systems with simple algebraic structures. First, we 
have applied the quasi-spin squeezed state to the pairing model by two approaches. One is 
dynamical approach, namely, the time-dependent variational approach. In this approach, 
the total fermion number is conserved automatically. We have imposed a consistency condi- 
tion for the rotation in the phase space. As a result, the ground state energy is realized up 
to 0(1/ Q 2 , 1/N 2 ), which is better than the result of the coherent state approach, namely, 
the Hartree-Fock-Bogoliubov approximation. The other is to obtain the ground state energy 
directly in the energy minimum condition with the chemical potential. We have estimated 
the ground state energy numerically in the variational method, and compared the result with 
that of the coherent state approach and the exact energy eigenvalue. It has been shown that 
this quasi-spin squeezed state approach outrages the coherent state approach in all regions 
of the strength of the paring interaction. 

Secondly, we have applied the quasi-spin squeezed state to the 0(4) model which has a 
su(2) x su(2) algebraic structure. Since we have the su(2) quasi-spin squeezed state to take 
account the pairing correlation, we have adopted the direct product of two su(2) quasi-spin 
squeezed states as a trial state. However, it has turned out that the quasi-spin squeezed state 
approach using the direct product of two quasi-spin squeezed states has not given good results 
compared with the 0(4)-coherent state approach. The reason is that the pairing correlation 
can not be included appropriately, and the direct product of two quasi-spin squeezed states 
is not reduced to the quasi-spin squeezed state for pairing model when \ — 0. Thus, we have 
constructed another squeezed state free from the 0(4) algebra, which we call the fermionic 
squeezed state. This improved squeezed state includes partially the quadrapole correlation 
and is reduced to the quasi-spin squeezed state when two parameters are identical. Although 
the fermionic squeezed state has given rather good results than the direct product state, it 
may be not enough to include the quadrapole correlation when the quadrapole interaction 
is dominated (G < 2x)- This is the further problem to find more suitable state. 

For further application of a quasi-spin squeezed state approach, it is interesting to inves- 
tigate a nuclear sii(2)-model which interact with the environment represented by a harmonic 
oscillator. 0* In this model, a certain case, a dissipative process has been realized. However, 
in our previous treatment, a quantum fluctuation has not been took into account. There- 
fore, a quasi-spin squeezed state approach might be suitable to introduce and investigate the 
effect of quantum fluctuation into the system. 
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Appendix A 

Various expectation values by using the direct product of two 

quasi-spin squeezed states 

In this appendix, we summarize the various expectation values with respect to the direct 
product of the quasi-spin squeezed states. We introduce new notations as ^(/3r/3|) = 
#{PoPo) = *o, #1 = dVi/dlfal 2 and % = d^ /d\(3 \ 2 . 



MjBA 



\M B \^ 




\BM) = 




WM) = 




WM = 


16|/3/| 


\M 2 \^) = 


16|/?o 


\B\Bj\^) 


= 2|/3/ 



B ] M \ 



2 ^+16|/3 |4 1 



^\B M \ 



\Bo\i/> 
b\M \ 



\B. 



\m i bM = o, 

, 



o 



1 \ % 4|/3/| 4 ^' 



BIB, 



BjB}\ 



BoBll 



2|/3o| 2 



25/7 ^/ 25/ >V 
1 \ *U 4|/? | 4 ^ 



o 

2S Q J &o 



O 



2S &o 



\B\BM)+ 1 



IM/I 



\bIb, 



o°o\ 



1 - 



25/ 

'MMo\ 



2S, 



o 



2/3 *'° 



\B 



Bh 2 \^)=2(3* ^ 1 



(A-l) 
(A-2) 

(A-3) 



(A-4) 



(A-5) 



(A-6) 



(A-7) 
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Thus, the expectation value of the Hamiltonian in ( |5-1|) can be calculated by using the above 
expectation values. 



Appendix B 

Various expectation values by using the fermionic squeezed state 

In this appendix, we summarize the various expectation values with respect to the 
fermionic squeezed state developed in §5-3. 
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n=0 



r-l 

*n(i 

p=i 



(2n-r-l) 
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9=1 



P 



2$ i J V 2 S () 



Q 



(s\M \s) = U E E l&l* IA| 4 - 2r 



2r 



r! (2n — r)! 



r-l 

xn(i 

P =i 



25 
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25/ 
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25 7 
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1 I n=0 r=0 
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25, 
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n i 
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2S T 
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The expectation values of other products of two operators, except for the complex conjugates 
of the above operators, are zero. 
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